Geometric phase of a two-level system in a dissipative environment 
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The geometric (Berry) phase of a two-level system in a dissipative environment is analyzed by 
using the second-quantized formulation, which provides a unified and gauge-invariant treatment of 
adiabatic and nonadiabatic phases and is thus applicable to a quantitative analysis of transitional 
regions away from ideal adiabaticity. In view of the recent experimental observation of the Berry 
phase in a superconducting qubit, we illustrate our formulation for a concrete adiabatic case in 
the Ohmic dissipation. The correction to the total phase together with the geometry-dependent 
dephasing time is given in a transparent way. The behavior of the geometric phase away from ideal 
adiabaticity is also analyzed in some detail. 
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I. INTRODUCTION 

A superconducting nanostructure [Ij] with its poten- 
tial scalability leads to a promising solid-state platform 
for quantum information processing [3, 01 ■ The coherent 
control of macroscopic quantum states in superconduct- 
ing circuits [1, in particular two- level quantum sys- 
tems, makes it possible to observe the geometric phases 
(GPs). If the evolution is adiabatic, CP is known as the 
Berry phase (BP) [^, which arises from the cyclic evolu- 
tion of a quantum mechanical system and depends only 
on the area spanned in the parameter space [7| . The ex- 
periments to observe BP have been implemented recently 
by using the superconducting qubit ^, iQj . It is intrigu- 
ing that a geometry-dependent dephasing was detected 
Q and it in general indicates a coupling of the qubit with 
its environment, which plays a non- negligible role in the 
superconducting circuits Given the argument that 

BP has an intrinsic fault-tolerant robustness and is po- 
tentially used as the geometric quantum logic gate 
the environmental effects on BP need to be quantitatively 
investigated. This issue of intrinsic robustness is very im- 
portant, but it appears that no consensus on this issue 
yet [11. 

The effects of environment on GPs have been analyzed 
by several authors from various points of view: for exam- 
ple, the classic field fluctuation [T3| , the quantum jump 
[l4| . GP distributions [1^, etc.. For macroscopic quan- 
tum states of superconducting circuits, the behavior of 
environment is generally simulated by an infinite num- 
ber of harmonic oscillators with some definite spectrum 
distribution of frequency [l^. Specifically, for a single 
qubit interacting with its environment, it is effectively 
described by a spin-boson model, which has been used 
in the past to analyze the quantum decoherence due to 
dissipation [l^ . This model can also be used for the 
analysis of GPs in dissipative environment and this di- 
rection has been initiated in [3 [H ■ 

In this paper we exploit a general method for treatment 



of the environmental effects in the interaction picture as 
a correction to effective energy eigenvalues appearing in 
the evolution operator. We here recall that the geometric 
phases, either adiabatic or nonadiabatic, are associated 
with the time development of the state vector typically 
during one cyclic evolution. The time development of the 
state vector is described by the Schrodinger equation, and 
the time development is entirely generated by the Hamil- 
tonian. The evolution operator, which is generically de- 
fined by \4'{t)) = J7(t, O)IV'(O)), if evaluated exactly thus 
contains all the information about the geometric phases. 
The primary object of our study is thus the evolution 
operator. 

To the second-order perturbation, we give a formula for 
studying the effects of the environment on GPs for a spin 
in the magnetic field at zero temperature, T ~ {). Com- 
pared with the previous relevant work 'l8], our formula- 
tion is valid for adiabatic [6] and nonadiabatic (cyclic or 
non-cyclic) [23, Hil evolutions. It shows the effects from 
a different point of view and is easily extended to analyze 
the dissipation-related issues such as adiabatic quantum 
computation (2^ and quantum decoherence [23| . In the 
following we first derive the general formula and then 
take the adiabatic limit to analyze the environmental ef- 
fects on BP in a concrete manner. The modification of 
GPs by environmental effects in the regions away from 
ideal adiabaticity is also analyzed. 

This article is organized as follows. In Sec. [TTl we 
review the second-quantized formulation of GPs in the 
absence of environment, which is exact for a single-spin 
system in magnetic field. Based on the formulation, we 
take the effects of environment into consideration in Sec. 
mil in which a perturbation correction to GPs originat- 
ing from the coupling with environment is displayed. As 
illustrations, in Sees. IIII Al and IIII Bl we show the effects 
of environment on adiabatic Berry phase and nonadia- 
batic geometric phase, respectively. Finally, Sec. IIVI is 
dedicated to conclusion. 
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II. REVIEW OF THE SECOND-QUANTIZED 
FORMULATION OF GEOMETRIC PHASES 

In this section, let us review the basic idea of the 
second-quantized formulation of GPs without taking dis- 
sipation into account. It was introduced by one of the 
present authors in [l^, HB, [H, US]- Instead of describ- 
ing a general theory, we consider the GPs of a spin in 
rotating magnetic field to demonstrate how the method 
works. Note that such a situation was already discussed 

in m,!!!. 

We denote a rotating background (magnetic) field by 
B{t) = B sin 9 cos (f {t), sin 9 sin If (t), cos 9) and (p{t) = 
loqI with a constant angular velocity loq. The action for 
the system in the second-quantized formulation (h = 1) 
is 



S = 



(1) 



where the field operator is expanded as 'ip{t,x) = 
c;(t)u;;(i) with the anti-commutation relation, 
{ci(<), cj„(t)} — Sim. For later use, here we define the 
Fock states by \l) — C;(0)|0) with the vacuum state |0) 
satisfying c;(0)|0) = 0. We use only the very elementary 
aspect of the second quantization of the fermion to clarify 
the hidden gauge symmetry that controls all GPs. 

For the above specific magnetic field with time- 
independent 9, the effective Hamiltonian for the isolated 
spin system is exactly diagonalized and, consequently, 
the isolated spin system is exactly solvable if one chooses 
the basis vectors as 



w+it) = 

with i9 = 
by 



-Mt) CO 
sm^ 



— 9o and the constant parameter 9q defined 



tan 9o 



ujQ sm ( 



B + ujo cos 9 



(3) 



Then we have 



wl.{t)hw±{t) — ^-B cos9o, 

wiit)idtw±it) = ia;o[l ± cos{9 ~ 9^)\ (4) 

with h ~ —B{t) ■ (t/2. In the operator formulation of the 
second-quantized theory, we obtain a diagonalized effec- 
tive Hamiltonian, Hcs{t) = J2i=± Eic\{t)ci{t), where two 
time-independent effective energy eigenvalues are given 
by 

E± = wi{t'){h-idt')w±{t') 

= T^Bcos9o-^LOo[l±cos{9-9o)]. (5) 
By noting the Heisenberg equation of motion 



it is confirmed that one can write 

ci{t) = uHt)ci{0)U{t) 

by introducing the "Schrodinger picture" effective Hamil- 
tonian 7icff(i) = Eicl{0)ci{0) and the second- 
quantized formula of the evolution operator defined by 



U{t) = rexp[- 



-H,s{t')dt'] 



(6) 



where T represents a time-ordering product. In general, 
we have Hcsit) = Z);.™ ^i,m(*)4(0)cm(0) for the time- 
dependent h{t), and the adiabatic approximation corre- 
sponds to an approximate diagonalization of Hcsit). See 
IM US m, [23 for more details. 

For the Schrodinger equation idtilj±{t) = h'ip±{t) with 
initial condition V'±(0) = ^"±(0), its exact solution is 
given in the second-quantized notation [2^, [l^] 

^±(<) = (0|^(i)4(0)|0) 



w± (t) exp 



'i / dt'wl{t'){h^idt')w±{t') 
Jo 



(7) 



where the exponent has been calculated in Eq. ([5]). Since 
w±{T) = w±(0) with the period T — 27r/|Li;o|j the solu- 
tion is cyclic [20] and, as an exact solution, it is applicable 
to the nonadiabatic case also. For an arbitrary time- 
dependent B(t), any exact solution of the Schrodinger 
equation can be written in the last form of Eq.([7|) , if 
one chooses basis vectors w±{t) suitably ^1]. But the 
periodicity w± (T) = w± (0) is generally lost and thus the 
solution becomes non-cyclic [21|. 

Actually, at the adiabatic limit IciJo/^l <C 1, 6*0 in 
Eq. ([3]) approaches zero so that the conventional BP 
7r(l ± cos9) [1] is recovered from Eqs. (j4|) and On 
the other hand, at the nonadiabatic limit |wo/-B| ^ 1, 6*0 
approaches 9 so that GP in Eq. ([7]) vanishes. Namely, 
the adiabatic BP is smoothly connected to the trivial 
phase inside the exact solution We can thus ana- 

lyze a transitional region from the adiabatic limit to the 
nonadiabatic region, which was not possible in the past 
formulation. 

One can assign a gauge-invariant meaning to the GP 
under general adiabatic or nonadiabatic evolution. To 
see this, let us recall that the field variable i[j(t,x) — 
X]/=± in Eq.([T]) is invariant under the simulta- 

neous replacements [25[ 

Q(t)^e-'"'(*)Q(i), wi{t)-.e'^'^'^wi{t), (8) 

and thus basic action ([T|) [even with dissipation; see 
Eg. pop ] is invariant under this exact gauge symmetry. 
One then confirms that the exact Schrodinger amplitude 
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tPi(t) = {0\'ip{t)cl {0)\0) in Eq.dT]) is transformed under 
this gauge symmetry as ~^ exp [jQ;;(0)]-!/'i(i) inde- 

pendently of t. The product ipj {0)ipi{t) is thus manifestly 
gauge-invariant. Its phase after subtracting the gauge- 
invariant "dynamical phase" (DP) dtw\ {t)hwi{t) be- 
comes 



A =argJ wJ(OVz(T)exp 



dtw\{t)idtWi {t) 



which is also manifestly gauge-invariant. This Pi is un- 
derstood as the holonomy of the basis vector associated 
with exact hidden local symmetry ([8]) for all GPs, ei- 
ther adiabatic or nonadiabatic, as explained in detail in 
[2^ . This construction is a generalization of BP for the 
generic case Ei ^ 0, for which the Schrodinger aniplitude 
does not satisfy the parallel transport condition [28] but 
the basis vector can satisfy the parallel transport condi- 
tion with the help of gauge symmetry ([5]) [26[. For the 
noncyclic case, one can still identify Eq.([9]) as a gauge- 
invariant noncyclic GP 21]. 

We here briefly compare the above unified formula- 
tion of GPs to the conventional formulation where the 
adiabatic phase is defined to be invariant under the 
symmetry identical to the above hidden symmetry ([8]), 
whereas the nonadiabatic phase is defined to be in- 
variant in the so-called projective Hilbcrt space with 
the equivalence class {e'"(*V(t)} [20|, |2l|. As a con- 
sequence, the gauge invariant nonadiabatic phase (3 — 
arg{^/>t(0)V'(T)exp[i/J'd#t(t)i5^^(i)]} is non- 

local and nonlinear in the Schrodinger amplitude i-!{t), 
which causes certain complications as was noted by Mar- 
zlin et al. [13]. In contrast, our j3i in Eq. ([9]), which nu- 
merically agrees with the Aharonov-Anandan /3 when one 
uses exact solution ([7]) in /3, is bilinear in the Schrodinger 
amplitude. 



III. GEOMETRIC PHASES IN DISSIPATION 

Coming back to the qubit in a noisy environment, ac- 
tion ID) after taking the environment into consideration 
is written in the second-quantized formulation as 

S = j dt['ij^[t)(i^^+B-cT/2)^{t)+Y,(pJc. 
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(10) 



where the environment is effectively described by an in- 
finite number of bosonic oscillations [s^ . Note that we 
work on the case of the vanishing temperature T = in 
the present paper. Accordingly, the effective Hamilto- 



nian after considering dissipation turns to be 



(11) 



7 (9) where Xa — 



.{da - dl), p. 



{da ~ di){wj^a^wi)cl^ci, 



and the time-independent effective energy eigenvalues Ei 
have been given in Eq. ([5]). The exact state vector for 
the qubit with the initial condition ijji (0) — wi (0) is given 
byV'/(i) = (0lV'(t)4(0)|0),i.e., 

Mt) - Vw™(t)(mirexp[-i / n,s{t')dt']\l) 

JO 



J2w^{t){m\U{t)\l). 



(12) 



In our specific example in Eq. (jlip . Ti.cs{t) is time in- 
dependent and simplifies calculations. After integrating 
out the environmental freedom, Eq. (jl2p assumes the 
form 



E 

m— lb 



m\ cxp < — ^ 



J2 i?fc4(0)4-(0) 



fc=± 



J2^kk'{t')cl{0)ck'{0) 



k,k' 



dt'\\l)w^{t), (13) 



in which the second term in the square brackets is the 
lowest-order "self-energy" correction, to be evaluated be- 
low due to the interaction with environment. 

Actually, our primary object of interest is evolution op- 
erator defined in terms of effective Hamiltonian pip . 
One can then integrate out the bosonic freedom, which 
is Gaussian, exactly in the path-integral representation 
of the evolution operator. 



Vxa'Dpa'Dip'Dip expliS"}, 



(14) 



where the action S is given by Eq. (jTO|) . This proce- 
dure is somewhat analogous to the Gaussian integral of 
the time component of the electromagnetic field Aq in 
quantum electrodynamics defined by the Coulomb gauge. 
One then obtains a four-fermion coupling analogous to 
the Coulomb interaction. But in the present case the 
Coulomb potential is replaced by the bosonic free prop- 
agator defined by (see Eq. (3) of [31*]) 

y2{TCaXa{t)C0Xp{t')) = - f ^ dw'J{uj') 

duj f i i \ ^_,^(t_t') 



27r V w — w' -f ie lo -\- lo' — ie 
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FIG. 1: The Feynman diagram for the self-energy correction 
"Smk in the one-loop order: the solid line stands for the spin 
and the wavy line stands for the bosonic freedom. 



where we replaced the summation over a in the x prop- 
agator by an effective spectral density, 



maOJo 



This spectral density typically has a power-law behav- 
ior at low frequencies ,17]. Of particular interest is 
the Ohmic dissipation, corresponding to a spectrum 
J(a;) = r]uj, which is linear at low frequencies up to some 
high-frequency cutoff (^c > B). The dimensionless 
parameter ry reflects the strength of dissipation. Here 
we concentrate on weak dissipation, yy <C 1, since only 
this regime is relevant for quantum-state engineering [j] . 
One may then perform the fermionic path integral corre- 
sponding to the lowest-order perturbative correction to 
the fermion self-energy depicted in Fig.l. Alternatively, 
one can perform the same calculation by using the Dyson 
formula S — Texp[— i J Hi{t)dt\ in the interaction pic- 
ture. 

The actual evaluation of the Feynman diagram in Fig. 
[T] for the Ohmic case is straightforward by following 
the second-quantized formulation of the Caldeira-Leggett 
model [31] . The self-energy correction in the one- loop or- 
der is then given by (see Eq. (5) of [31]) 



-(1) 



vJ2^wla,wi){w}a,Wk) [i{E - Ei)Q{E - Ei 
i=± 
{El 



E) 



In 



E-Ei 



(16) 



in which O is the step function. The first term in Eq. 
(I16p . which is imaginary, gives the decay width of the 
fcth level as 



^ i=± 



wla,wi\\Ek^Ei)eiEk-Ei) 



(17) 



which vanishes for the ground state due to the step func- 
tion. It indicates that the excited state decays to the 
ground state by emitting soft bosonic excitations. Below 
we will see that it characterizes the dephasing time scale 
of the spin system in magnetic field. The second term 
of Eq. (|16p is a correction to the effective energy. As 
a result, the total effective energy, to the order of 0{r]), 
becomes 

V 



^tot 



Ek + -y^\wla,wi\ {Ei~Ek) 



X In 



i=± 



Ek - El 



(18) 



The Schrodinger amplitude in Eq. (|T3|) to the accuracy 
of the lowest-order correction 0(77) is thus given by using 
Eqs.HH) and ((T8|) as 



Ef°\t')dt']wi{t). (19) 



It appears that the probability conservation for the 
higher-energy state is violated in Eq. (fT9|) . Mathemati- 
cally, this arises from the fact that we evaluated the "per- 
sistent" amplitude for the single-spin state under the in- 
fluence of dissipation. Evolution operator ^ is formally 
unitary in the present case also. Thus the probability 
should be preserved in our formulation and also in the 
formulation with the density matrix in the Appendix. 

In fact, the unitarity of the evolution operator, 
U'^{t)U{t) = 1, is preserved if one evaluates 

wi^{t){l+\U{t)\l+) + wi_{t)^B{t){l.;hosoiis\U{t)\l+), 

(20) 

since 'Wi_{t){l^\U{t)\l+) = in the present model, where 
{l±\ respectively stand for the higher- and lower-energy 
states of the spin, and (/_; bosons] stand for the final 
states of the lower energy spin state together with soft 
bosonic excitations. (To be more precise, one may write 
{l±\ as {l± \ ® (0| and (L; bosons] as (/_| ® (bosons] with 
(0| standing for the bosonic vacuum. But we use the sim- 
plified notation in this paper.) Note that the description 
in terms of effective spectral density (jlSp is different from 
the description in terms of a definite number of bosonic 
quanta. In the spirit of the Caldeira-Leggett model p^ . 
we do not assign a physical significance to the soft bosonic 
excitations [SQj. In fact, the actual cause of the dissipa- 
tion could be completely different from an ensemble of an 
infinite number of harmonic oscillators, though we con- 
sider that those harmonic oscillators, if suitably chosen, 
can mimic the actual dissipation. The presence of the 
dissipation is thus manifested by the decrease in the per- 
sistent amplitude w;^ (t)(Z+|C/(t)|Z+). [In this respect, our 
bosonic freedom is different from the time-component of 
the electromagnetic field in QED, which does not infiu- 
ence the unitarity. In path integral (fTi|) . one needs to 
add a Schwinger's source term to the bosonic freedom 
also to analyze the unitarity.] 

Since the final states wi_{t)ipB{t){l~\hosous\ are or- 
thogonal to both of w;^(t)(Z+| and 'Wi_{t){l-\ in Eq. 
(pn| and thus do not interfere with these states when 
one considers only the observables of spin freedom as 
in the present case, the decrease in the persistent am- 
plitude wi^{t){lj^\U{t)\l-[.) is also understood as an in- 
dicator of the quantum decoherence. See Eq. (|A3p in 
the Appendix. The inverse of the decay width in Eq. 
(ITTt is thus related to the two characteristic time scales 
of the qubits, namely, T^ciax and Tip [l[. The relaxation 
time scale Troiax is the time required for a qubit to relax 
from the first excited state to the ground state, involv- 
ing energy loss. The dephasing time scale is the time 
over which the off-diagonal [in the preferred basis vectors 
w±{t)\ elements of the qubit's reduced density matrix de- 
cay to zero in the formulation with the density matrix as 
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FIG. 2: (Color online) Berry phase 4>J^^ for various values of 
cutoff uic/B. The black line corresponds to the dissipation- 
less case, while marker points represent the case with dissi- 
pation. We set the dimensionless parameter r] = 0.3 and take 
sign(tt;o) = -1-1. The values of the period T — 2tt/uio and the 
field strength B/2ir = 50 x 10''(Hz) are taken from [1]. 



Here we mention that to derive above t~ , we first set 
^0 — sin 9loq/B and then let loq/B — > in formula 
([T7|) . For \ujo\ <C -B < oJc, the logarithm in the last 
term is nearly independent of geometry and its sign is 
very important because it determines whether the cor- 
rection term is added or subtracted. After evolving one 
cycle T (27r/_B ^ T ^ t^), if one defines the phase 
= $pP -I- $BP ^ jT ^tot(^)^^ 27rn (n =integer), 

we may identify DP and BP respectively as 



^-BT± -sin^ OBTln 



1 ± 



B 



sgn(ijjo) ±r2 =F 27]sin^ (?cos(?ln 



1± — 

B 



where the solid angle Q. — 7r(l — cos 61). The first terms 
in both DP and BP above are the ones in the absence 
of noise, while the second terms arise from the coupling 
between qubit and environment. Our results on the cor- 
rection of BP and the geometric dephasing factor are 
consistent with those found by Whitney and co-workers 

MM- 

If LOc < 2B, the BP correction of the ground state 
would have the same sign as that of the excited state. 
Otherwise, the sign will be opposite. For an illustration, 
we plot <I>BP with respect to the polar angle 6 for various 
values of the cutoff uUc/B in Fig. [21 

In passing, for a super-Ohmic case J{uj) = ijuj^ the 
dephasing time scale is given by r"^ — rjsir? 9B'^{B + 
loqCOsO) and the correction to BP becomes A$bp — 
Ts\gn{ujo)'2ri sin^ 9 cos e{uoll 2 - B'^\n\l±uOc/B\) in the 
adiabatic limit. 



in the Appendix. In the present model, these two time 
scales are of the same order of magnitude, and we choose 
= 2/r/+ by using the decay width in Eq. (fT7)) . We 
thus need to satisfy T <C or 



(21) 



to make a sensible measurement of geometric phases, 
which are defined for pure states. The gauge-invariant 
geometric phases for mixed states can be defined, but 
their measurement is generally rather indirect (26j . 



A. Adiabatic Berry phase 

We now come to the analysis of GPs with dissipation. 
For dcfiniteness, let us first analyze the environmental 
effects on the adiabatic BP. In the adiabatic limit 9q ~ 
sh\6L0f)/B — !■ in Eq. ([5]), one has the non-vanishing 
quantities 



rj s\T? 6{B — ojQ cos 9), 

t]-B -\uoo{l±cos9)±- sn? 

2 2 TT 



-Wo cos ( 



I In 



1 ± 



B + loq cos I 



B. Nonadiabatic geometrical phase 



We emphasize that our formula p9)) is valid for both of 
the adiabatic and nonadiabatic cases. Thus the modifica- 
tion of nonadiabatic (Aharonov-Anandan) GP due to the 
dissipation is analyzed quantitatively. A technical issue 
involved here is that a clear separation of dissipation- 
induced GP from dissipation-induced DP is not possible; 
One already recognizes this tendency even in the case 
without dissipation in Eq. We here tentatively em- 
ploy the following procedure for identifying the geometric 
phase away from the ideal adiabaticity with dissipation. 
As for the phase without dissipation, we adopt the sepa- 
ration in Eq. (|4]). For the phase induced by dissipation, 
we define the geometric part by 

TA4"*) = T[Ai?fe(L^o) - l^Ekiuoo = 0)], 

where /S.Ek{LOo) stands for the second term in Eq. (fTS]) 
proportional to 77, and T = 27r/|ajo|- This identification 
may be reasonable if both of loq and 77 are small. Thus 
the total geometric phase is 
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Jl/4 Jl/2 3ll/4 Jl 

e 



FIG. 3: (Color online) Non-adiabatic geometric phases for dif- 
ferent evolution rates uq/B, in the presence or absence of envi- 
ronment. All lines including solid and dashed ones correspond 
to the nonadiabatic GPs without environment coupling (i.e., 
flk), while all marker points represent GPs, j'£;(g'=°™<=') ^ with 
environment coupling. The adiabatic limit of |a;o/-B| <C 1 
is shown by black color. We set the dimensionless parameter 
rj = 0.3 and take sign(tiJo) = +1- Other parameters are chosen 
as BI2-K = 50 X 10^ (Hz) and Uc = ZB. 



tends to be restricted in the neighborhood of zero, the 
trivial value of GP. 

In the transitional region from adiabatic to nonadia- 
batic limit, the effect of environment on GP changes as 
increases and keeps the same pattern for diverse ujq/ 
as indicated by the red and blue lines in Fig. [31 The red 
dashed lines, corresponding to a relatively small value of 
lhq/B, tends to recover the adiabatic case shown in Fig. 
[21 while the blue dashed lines away from both the adia- 
batic and nonadiabatic limits reflect the transitional re- 
gion. One may note that at 6* = tt GP takes one of values 
of 0, ±7r, and ±27r, which sharply depends on the value of 
ujq/B. This can be easily understood from the definition 
of ^0 in Eq. ^ . In addition, it should be mentioned that 
for a fixed value of cvo/B the Ti;|^™'°°*^ and TE^^'°'^°'^ 
are not symmetric about the axis of T^^^^scomot) _ q 
to the existence of the logarithmic term in Eq. (|18p . As 
for the role of the coupling strength 77, it determines the 
magnitude of deviations, while having little effect on the 
pattern of deviations at least in the perturbative domain. 

Finally, let us briefly analyze the geometric dephasing 
factor for the nonadiabatic case. From Eq. p?)) . the 
nonvanishing decay width is given by 

-r_ =T]sm^{B-BQ)[B cos 00 + ^0 cos{6 — Oq)]. 

Its geometric dependence is complicated because of the 
involvement of 6*0 through Eq. ([3]). At the adiabatic 
limit, it returns to the discussed in Sec. IIII Al At 
the nonadiabatic limit, we have {0 — Oq) ~ and thus re- 
alize a dephasing-free situation. Unfortunately, the limit 
corresponds to the trivial geometric phase. 



which is illustrated in Fig. [31 with respect to the vari- 
able 9. Here ftk — ^^Bq stands for the nonadiabatic 
GP without dissipation obtained by multiplying the pe- 
riod T to the second "geometric energy" in Eq. ^ and 
rigj, = n[l — cos{9 — 9q)] stands for the (nonadiabatic) 
solid angle with 9q defined in Eq. ([3]). When taking the 

adiabatic limit of Iti^o/^l *C Ij T£'^^''°™'^*'' approaches 
the adiabatic Berry phase with dissipation. When tak- 
ing rj ~ 0, ig reduced to the nonadiabatic GP 
rik = if^So without dissipation. 

Figure [3l illustrates how the observed GPs, with and 
without the effects of environment, deviate from the sup- 
posedly topological BP, ±7r(l — COS0), when the variables 
ujQ and T] deviate away from negligibly small values. It 
can be seen that the nonadiabatic GP is affected sensi- 
tively by the change in the dynamical parameter loq/B, 
while the adiabatic BP is supposed to be immune to small 
fluctuations of ujq/B. This is related to the basic issue of 
the intrinsic robustness of BP. When ojq/B becomes big 
{e.g., > 1 in our choice of the parameter values), the non- 
monotonous magenta dashed lines localized near zero in 
Fig. [21 show that the dissipationless nonadiabatic GP fifc 



IV. CONCLUSION 

We have presented a gauge-invariant formulation of 
GPs for a two-level system in dissipative environment, 
which is applicable to both of adiabatic and nonadia- 
batic cases. Our formulation may be useful for under- 
standing the experimental observation of BP such as in a 
recent superconducting qubit [8] and will provide a start- 
ing point for the future quantitative analysis of the basic 
issues of intrinsic robustness and the behavior of GP 
in the transitional region away from ideal adiabaticity. 
The analysis of the region away from ideal adiabaticity 
is expected to be crucial in any practical application of 
geometric phases. 
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APPENDIX A: DENSITY MATRIX 

One may be interested in the mixed states which are 
described by a density matrix in general. The time de- 
velopment of the density matrix p is described by the 
evolution operator as 

p{t)^Uit)p{0)U\t), (Al) 

at the vanishing temperature T = 0. The evaluation of 
U{t) is thus sufficient. For our purpose of simulating the 
dissipation by an infinite number of harmonic oscillators, 
one may trace out p{t) = U{t)p{0)U'' {t) with respect to 
all those final states which contain bosonic excitations. 
In our simple model, the bosonic excitations are included 
only in the lower-energy state of the fermion as in Eq. 
(1^0)) . If one defines a state 

iP{t) ^aibi^{t)+b^i_{t), (A2) 

with constants a and b, which is no more cyclic even 
without dissipation since ip{T) ^ ?A(0) up to a phase, 
and starts with a pure state p{0) — |i/'(0))(V'(0)|, one 
obtains the mixed state after partial tracing of p{t) over 
those final bosonic excitations as 

m - \m){m\ + w\'[HiAo)\'-\^iAt)f]^ 

\^ijt)){i^i^{t)\, (A3) 



with |'!/'i+(0)P = 1. We here used the unitarity relation 
inEq. jfOl), 



bosons 

= l^;.(0)P, (A4) 

with = 1. The first term in Eq. (|A3|) arises 

from the trivial bosonic vacuum and the second term in 
Eq. (IA3p from the nontrivial bosonic states. The density 
matrix p{t) approaches 

m^{W\' + \b\'Mi-{t)){i:^_{t)l (A5) 

for t large. Thus the total trace Tip{t) = |ap -1- |&p 
is preserved during the time development in the present 
normalization of the density matrix. 

One may define gauge invariant geometric phases for 
the mixed state described by this density matrix p{t) 
following the general formulation in (26j . An interesting 
case is obtained if one sets 6 = and considers the time 
interval to < t < to + T with to ^ for p{t). 



[1] J. Clarke and F. K. Wilhelm, Nature (London) 453, 1031 
(2008). 

[2] J. Majer et al.. Nature (London) 449, 443 (2007). 

[3] M. A. Sillanpaa, J. I. Park, and R. W. Simmonds, Nature 

(London) 449, 438 (2007), and references therein. 
[4] Y. Makhlin, G. Schon, and A. Shnirman, Rev. Mod. 

Phys. 73, 357 (2001). 
[5] D. Vion et al.. Science 296, 886 (2002). 
[6] M. V. Berry, Proc. R. Soc. London, Ser. A 392, 45 (1984). 
[7] J. Anandan, Nature (London) 360, 307 (1992). 
[8] P. J. Leek, et al.. Science 318, 1889 (2007). 
[9] M. Mottonen, J. J. Vartiainen, and J. P. Pekola, Phys. 
Rev. Lett. 100, 177201 (2008). 
[10] Y. Nakamura, Yu. A. Pashkin, and J. S. Tsai, Nature 

(London) 398, 786 (1999). 
[11] P. Zanardi and M. Rasetti, Phys. Lett. A 264, 94 (1999). 
[12] A. Blais and A. M. S. Tremblay, Phys. Rev. A 67, 012308 
(2003). 

[13] G. De Chiara and G. M. Palma, Phys. Rev. Lett. 91, 
090404 (2003). 

[14] A. Carollo, I. Fuentes-Guridi, M. F. Santos, and V. Ve- 

dral, Phys. Rev. Lett. 90, 160402 (2003). 
[15] K.-P. Marzlin, S. Ghose, and B. C. Sanders, Phys. Rev. 

Lett. 93, 260402 (2004). 
[16] A. O. Caldeira and A. J. Leggett, Phys. Rev. Lett. 46, 

211 (1981). 

[17] A. J. Leggett, et al.. Rev. Mod. Phys. 59, 1 (1987). 



[18] R. S. Whitney and Y. Gefen, Phys. Rev. Lett. 90, 190402 
(2003). 

[19] R. S. Whitney, Y. Makhlin, A. Shnirman, and Y. Gefen, 

Phys. Rev. Lett. 94, 070407 (2005). 
[20] Y. Aharonov and J. Anandan, Phys. Rev. Lett. 58, 1593 

(1987) . 

[21] J. Samuel and R. Bhandari, Phys. Rev. Lett. 60, 2339 

(1988) . 

[22] E. Farhi et al.. Science 292, 472 (2001). 
[23] F. C. Lombardo and P. I. Villar, Phys. Rev. A 74, 042311 
(2006). 

[24] S. Deguchi and K. Fujikawa, Phys. Rev. A 72, 012111 
(2005). 

[25] K. Fujikawa, Phys. Rev. D 72, 025009 (2005). 

[26] K. Fujikawa, Ann. Phys. (N.Y.) 322, 1500 (2007). 

[27] K. Fujikawa, Phys. Rev. D 77, 045006 (2008). 

[28] B. Simon, Phys. Rev. Lett. 51, 2167 (1983). 

[29] K. Singh, D. M. Tong, K. Basu, J. L. Chen, and J. F. 
Du, Phys. Rev. A 67, 032106 (2003). 

[30] In the spirit of the Caldeira-Leggett model p^ . we treat 
the bosonic freedom as an auxiliary device to simulate 
the dissipation. We do not take the infinite number of 
bosonic excitations as physical objects. See K. Fujikawa, 
Phys. Rev. E 57, 5023 (1998). 

[31] K. Fujikawa, S. Iso, M. Sasaki, and H. Suzuki, Phys. Rev. 
Lett. 68, 1093 (1992). 



